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We study the process of multimode Spontaneous Parametric Down— Conversion (SPDC) in the lossy, one di- 
mensional waveguide. We propose a description using first order Correlation Functions (CF) in the fluorescence 
fields, as a very fruitful and easy approach providing us with a complete information about the final multimode 
state. We formulate the equation of the evolution of the multimode CF along the crystal using four character- 
istic length scales. We solve it analytically in the one mode case and numerically in the multimode case. We 
capture simultaneous effects of three wave mixing with ultrashort pump, linear propagation and attenuation, 
and we are able to divide the evolution into three stages and predict it qualitatively. We find that losses do 
not destroy the quantum properties of SPDC but stabilize the final state. 

Keywords: Multimode Spontaneous Parametric Down-Conversion; Squeezed states; Gaussian states of light; 
Dissipative dynamics; First order correlation functions. 

1. Introduction 

Optical Parametric Amplifiers (OPAs) are commonly used as a source of non-classic light. Oper- 
ated below threshold they produce squeezed light which can be used in a vast variety of quantum 
protocols, including state teleportation (0), quantum cryptography |J), interferometry (0) and 
precise measurements Q). 

A typical OPA comprises pumped nonlinear crystal placed in an optical cavity resonant at 
exactly half the pump frequency. In such systems narrowband, continuous-wave squeezing is 
obtained Q). In this paper we focus on the other possibility, when the nonlinear medium is 
pumped with a short pulse and a broadband squeezed light is produced in a single pass through 
the crystal in a process commonly called Spontaneous Parametric Down-Conversion (SPDC). 
This is especially efficient when a waveguide made in the nonlinear crystal is used to carry the 
pump and squeezed light (0). 

It is well known that optical losses destroy squeezing. However nonclassical character of those 
states still manifests itself in residual squeezing present even after severe attenuation. This be- 
havior is rather uncommon among nonclassical states of light. Most of them die rapidly subject 
to deleterious contact with environment. Therefore we have chosen to study the effects of dis- 
tributed losses, unavoidable in realistic nonlinear materials, on generation of squeezed light. It 
turns out that squeezing can still be generated in lossy media and it is even stabilized under 
certain conditions. However, the output state no longer saturates the uncertainty principle. 

In previous works the multimode state of squeezed light was described either with help of 
Green functions Q) or perturbative expansions Q) have been used. However those approaches 
cannot be used for numerical description of the lossy media, because additional, very large Green 
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functions appear. Therefore we have adopted a formalism which traces the evolution of first order 
correlation functions of the fluorescence fields. This approach is computationally effective and 
provides complete information about the final multimode state. Moreover the evolution of the 
correlation functions is simple and qualitatively easily predictable. Below we will describe various 
possible cases and assess characteristic widths of the correlation functions by pump pulse length, 
crystal thickness and its fundamental dispersive properties. 

An SPDC process in an ideally transparent nonlinear medium results in creation of multiple 
independent squeezed states (0, Q) . They occupy orthogonal modes of light field and can be in 
principle separated (0). Moreover each of the states is pure and saturates Heisenberg inequality 
for variance of x and p quadratures. 

However, as soon as loses come into play full description of the final state becomes much 
more complicated. The state can be cast into Williamson form (Til), which corresponds to a 
reconstruction as a multimode squeezed state subject to losses mixing components with various 
squeezing. Unfortunately there is no way to separate the field into uncorrelated modes. One 
can find a set of polychromatic quadratures of light Xi that are not correlated with any other 
commuting quadratures, however together with they conjugates they exceed minimal uncertainty 
bound. 

We describe the multimode state of light in terms of its experimentally important character- 
istics. We trace the evolution of the average number of photons, maximum squeezing and the 
product of variances of maximally squeezed quadrature and its conjugate. We differentiate three 
general stages of the evolution: an exponential growth and a squeezing, a linear temporal expan- 
sion without squeezing and a stabilization. We show that the boundaries between those stages 
typically correspond to characteristic lengths of either attenuation or pump-SPDC overlap. 

The paper is organized as follows. In Sec. 2. we introduce the correlation function and describe 
the adopted model of OPA. In Sec. 3. we present results of numerical simulations of the evolution 
of the correlation functions. In Sec. 4. we connect the correlation functions with the results of the 
Homodyne Detection (HD). By performing an optimization procedure we find the most squeezed 
mode of the light. Furthermore, we gather the results of our simulations and divide the evolution 
into stages. Finally Sec. 5. concludes the paper. 

2. Model 

In this section we develop a model of a lossy OPA. In particular we introduce correlation functions 
and equations for their evolution. We define four characteristic length scales characterizing the 
evolution of the SPDC light. 

We will focus on interactions in a waveguide made in a nonlinear medium Q) in which fields 
propagate along z axis in fixed spatial modes. We assume that the interaction of the pump at 
frequency 2cuq and the fluorescence at frequency cuq is perfectly phase matched. Due to the small 
intensity of light in SPDC the pump remains undepleted. We will first discuss a situation in 
which only the SPDC undergoes linear losses. Then we show how the evolution is modified when 
the pump is also damped. 

We begin by discussing the propagation of the field at a degereate frequency u>o in a lossless 
OPA. The equation of evolution of the SPDC annihilation operator a can be written as da/dz = 
a) /Lnl, (fill ) where characteristic length of the nonlinear coupling Ljvl is inversely proportional 
to the pump amplitude Pq and nonlinear coupling. Let us note, that this equation is formally 
identical to the equation of evolution of the classical field amplitude along the nonlinear medium. 

For a lossy nonlinear medium the description becomes more involved. Classical equation of 
motion would in such case contain a term —a/ La with La being attenuation length. In quantum 
case this term is inevitably accompanied by an additional vacuum noise: 
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Figure 1. Schematic representation of the model of lossy OPA. The OPA is divided into thin slices of thickness Az 

where f(z) are Langevin noise operators. They form a set of independent bosonic operators 
indexed by z, i.e.: [f(z), f* (z f )] = 5(z — z') representing vacuum modes (f(z)) = 0, (f^(z)f(z)) = 
(f(z)f(z)) = 0. The above evolution can be understood as an alternate application of squeezing 
and beamsplitter type losses, as illustrated in Fig. [TJ Eq. (|2.1h can be integrated to the following 
input-output relations: 



a(L) = e~ L/LA (a(0) cosh L/L NL + a t (0)sinh L/L h 

+ dr e ~ r/LA (f(L ~ r) cosh r/L NL + p(L - r) sinh r/L NL ) (2.2) 



Note that the solution involves entire set of f(z), although one expects that the final state on 
the output could be described without tracing the nuances of this dependence. Indeed, we may 
use the fact, that the state of SPDC is always Gaussian centered around the origin of the phase 
space. This is because we start from vacuum states that also have this property and the equations 
of motion (|2.ip are linear, corresponding to linear distortion of the phase space. Therefore the 
final state is completely described by its second order moments (d^d) and (ad). Both of them 
are c-numbers and their values can be easily traced by solving their evolution equations. They 
are in turn found by inserting Eq. (|2.ip into the z-derivatives of the second order moments: 
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Note, that in the above equations the Langevin noise operators f(z) do not appear, because 
their correlation functions are zero. The above equations can be solved analytically, completely 
describing evolution of the state of SPDC along the amplifier: 
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In principle the same result can be obtained from the Heisenberg transformation (|2.2j) . however 
this way we would first laboriously find the entire history of noise contribution only to trace it 
away in the next step. 

Let us briefly discuss this result, as it also captures important features of a full multimode 
evolution. For L A > Lnl/2 both (d^d) and (dd) grow exponentially. However, for losses exceeding 
critical level La < Lnl/2, sinh and cosh functions are dwarfed by the attenuation exponent and 
the (d^d) and (dd) ultimately reach steady state values. 
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When the pump is a short pulse and we wish to include material properties of the medium, 
a multimode description of the SPDC by a family of operators a(uj) becomes necessary. The 
dispersion of the crystal is reflected in the dependence of the wavenumber on frequency k(u). It 
is included in the propagation equations by adding a term of type ik(u)a(uj). Since the spectrum 
of SPDC in a waveguide is restricted by moderate phasematching bandwidth, we can expand 
k(uj) around SPDC central frequency ujq up to second order. The expansion coefficients are 
inverse phase velocity, inverse group velocity and second order dispersion respectively. Similar 
expansion around the pump frequency 2u)q yields its propagation parameters. We assume that 
the pump bandwidth is sufficiently narrow so that the effects of the second order dispersion on 
the pump pulse lengthening can be neglected. Due to phase matching assumption phase velocities 
of the SPDC and pump match and drop out. Therefore, to describe the material properties, it is 
sufficient that we know difference of inverse group velocities between the pump and the SPDC 
A/3i and dispersion coefficient for SPDC fa- 



Aft = — 
da; 



dk 
duj 



02 



d 2 k 



^=2^ duj2 



(2.5) 



We assume the spectral amplitude P{uj) of the pump pulse to be Gaussian corresponding to 
a duration rp and a maximum amplitude Pq 

PH = P e- (u - 2uiW/2 . (2.6) 

Now we can introduce a length scale Lov ° ver which pump pulse and the SPDC light overlap 
Lov = Tp/Af3i. Additionally one can define a crystal length Lr> over which the spectrum of the 
SPDC would get as narrowband as the pump spectrum, Ljj = rp/ (32 EL With those definitions 
at hand we can write a multimode equation of the propagation of a(co) in the reference frame of 
the pump pulse: 



da(uj) .( t p uo t p uj 2 \ I f , t , , 1 



The first term represents linear dispersive evolution in the pump reference frame. The second 
term describes nonlinear interaction, which couples SPDC components of frequency uj and uj' 
with pump components of frequency u + uj' . The last two terms describe losses. 

Since the above equation is linear its solution can be always expressed as an integral transform: 



a(u>, z) = I C(u, u')a(u', 0) dui' + J S{uj,oj')a) (uj' ,0) duj' 

+ j duj'dz'Cf(uj,uj',z')f(uj',z')+ / dio'dz' Sf(uj,uj' , z')p (uj' , z') (21 



The Green functions C(uj,uj'),S(uj,uj'),Cf(uj,uj',z'),Sf(uj,uj',z') contain all the information 
about the amplifier, in particular they are sufficient to describe the final state. However, finding 
numerical approximations of Cf(u),u)' ,z') and Sf(uJ,uj',z') would be very tedious and unneces- 
sary. The quantum state of SPDC is fully described by a multimode generalizations of the second 
order moments which become Correlation Functions (CF) (a^(uj)a(uj')}, (a(uj)a(uj')}. They ex- 
press the amplitude and the phase correlation of the light between two arbitrary frequencies uj 
and uj' . The evolution of the CF can be traced by solving their evolution equations. These are 



^ee Sec. 3 for detailed explanation of the role of Lp 
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found by differentiating (a* (uo)a(uo')) and (a(u)a(oo')) over z and using Eq. (|2 

g(at( ^° (U;/)) = (L» + L(a/)) <at(u;)a(a/)) - -?-<ot( W )a(a/)) 

+ — y — f doo" [P*(oo")(a(u - oo")a(oo')) + P(oo")(a(oo)a(oo" - J))*] (2.9a) 

PqLnl J 

d (a(U f ^ = {m + L(u/)) (o(oo)a(J)) - -L(o( W )a(o/)> + p-^-iV + 



+ — !— / <L/' 

PqLnl 



P(oo")(a)(oo" - u)a(u')) + P(oo")(a) (oo)a(uo" - to'))*] . (2.9b) 



Above L(oo) denotes the linear propagation operator in frequency domain L(ui) = irpoo/Lov + 
ir 2 oj 2 /(2L D ). 

Actually it is more intuitive to trace the evolution of the CF in the time domain, rather than in 
frequencies. It is meaningful to apply double Fourier transform to (a' (00)0,(00')) and (a(oo)a(oo')) . 
This way we obtain CF in time: 

(a\t)a(t')) = J doodJe^-^' 1 ' (a)(oo)a(oo')), {a(t)a(t')) = j ' doo doo' ' e -™t-m't> \a(oo)a(oo')) . 

(2.10) 



This transform can be applied to both sides of (|2.9p and one can obtain evolution equations for 
CF in time domain: 

d{aKt)a(t')) (r * , r s r U+ V (^\ < at (*) a (0> + P*(t)(a(t)a(t')) + m)a(t'))*P(t') 
= (A + Ctf) (aHt)a(t )) - ^ + _ 

(2.11a) 

a(a(y(Q) = 

oz La 

| P(t)(a\t)a(t')) + (aHt)a(t')yP(t') | P(t)5(t - Q 
PqLnl PqLnl 

The first term in (|2.1ip represents time shift and dispersive broadening, the second term repre- 
sents attenuation and the last term represents nonlinear interaction. The operators Ct represents 
convolution with a Fourier transform of L(oo) and can be rewritten as: 

r *^ r 2tp d , 2ir 2 P d 2 



L OV d(t + t') L D 8(t + t')8(t - t>) 

_ 2r P 8 it 2 ( d 2 d 2 \ 

1 + *' " L OV d(t + t>) L D \d(t + t') 2 + d(t - t') 2 J { ' 

where we on purpose changed the variables, so that differentiation along the diagonal and an- 
tidiagonal is emphasized. 



3. Evolution of the correlation functions 



Correlation Functions (CF) (af (oo)a(oo')) and (0,(00)0,(00')) completely describe the Gaussian quan- 
tum state of the SPDC. Therefore it is beneficial to understand evolution of the CF along the 
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crystal. Despite complicated form of the Eqs. (12. 9h . we will show that this evolution can be intu- 
itively understood and the CF shape is directly connected with four characteristic length scales: 
pump-SPDC overlap length Lov-, dispersive length Ld, attenuation length La and nonlinear 
length -Ltvl- Lov sets the effective length of the crystal over which the nonlinearity can build up 
SPDC in an undisturbed way. Ld is the characteristic length over which correlations between 
different moments in time arise. Finally La sets the crystal length over which contributions from 
the front slices of the crystal are erased and the state stabilizes. 

We have calculated the CF in the subsequent distances of the propagation by numerically 
solving matrix approximations of Eqs. (12. 9ft using the Split Step Fourier Method. We have verified 
that our simulation is accurate by repeating it on grids finer in each dimension. Calculations 
were carried out for a comprehensive set of the characteristic lengths in order to understand a 
mechanisms governing the evolution of CF. Our findings are most conveniently described in the 
time domain. 

Firstly, let us discuss the lossless case, i.e. La — ► oo. In the Fig. [2] we present results of the 
numerical simulations performed with the following parameters: pump-SPDC overlap length 
Lov = Lnl/W and dispersive length Ld = 3L^vl- Such choice gives an illustrative insight into 
a typical evolution of the CF inside the crystal when both a linear propagation and nonlinear 
coupling play comparable role. Using contour plots in the Fig. [2] we drew snapshots of the 
absolute value of the CF in the time domain \(a^ (t)a(t'))\ and \(a(t)a(t')}\. Snapshots were taken 
at crystal lengths L = Lnl/20, £jvl/3 and Ljvl- In the Fig. [3]we draw contours of the phase of 
the CF at L = Lnl- 

Let us analyze the plots presented in the Fig. [2j In the absence of losses the evolution is deter- 
mined by the interplay between two terms of Eqs. f|2. 1 1 [) : linear, responsible for the shifting and 
spreading of the CF and nonlinear, responsible for the growth of the CF. On a two-dimensional 
plane of t and t' the nonlinear terms act only in the region where P(t) and P(t') are nonzero. 
In particular the only seed term P(t)5(t — t') is located on the diagonal. Over distances when 
dispersion is negligible L <C Lov,Ld, the correlation functions grow on the diagonal accord- 
ing to single mode solutions (|2.4p . Then they diverge due to material dispersion. Hence, after 
L = Lnl/20 of the propagation (a(t)a(t')} arises around the diagonal pump seed region (the di- 



agonal of the square in the plots) as seen in the Fig. 2(b) At the same time the other correlation 
function (at (t)a(t')) is seeded by the term involving the pump and (a(t)a(t')). 

Let us note a finite antidiagonal width of both (at(t)a(t')) and (a(t)a(t')) . It is an effect of 
the second order dispersion of the crystal i.e. finite Ld. Indeed along the antidiagonal alone 
t + t' = const the linear part of the equation of evolution of (a(t)a(t')) (12. lib)) has a form of a 
diffusion equation d(a(t)a(t')) / dz oc d 2 (a(t)a(t')} / d(t — t') 2 with a diffusion coefficient inversely 
proportional to Ld- Therefore the (a(i)a(i')) spreads in the antidiagonal direction and its width 
can be estimated as t v \J 'Lj 'Ld- Within the pumped region where P(t) and P(t') are nonzero, 
this spread is then imprinted onto (at (t)a(t')) due to nonlinear term present in Eq. (|2,llap . Note 
that the antidiagonal width of the (at (t)a(t')) is inversely proportional to the diagonal width 
of (at (uj)a(uj')) . The latter is also the spectral bandwidth of SPDC which is now estimated as 
Tp 1 \JL~dJL for short crystal lengths. 



In the next two plots, Fig. 2(c) and 2(d) we capture CF after an evolution distance Lnl/3. 
The influence of the group velocity difference between the pump and SPDC becomes apparent. 
In Eqs. (|2.1ip it causes a shift of (at(t)a(t')) and (a(t)a(t')) along the diagonal with a speed 
of Tp/Lov- Both CF unwind from the pump seed region. Let us note that (a(t)a(t')) is under 
continuous influence of dispersion caused diffusion in the antidiagonal direction. It assumes 
triangular shape with broad part corresponding to photon pairs born in the front slices of the 
crystal and narrow part corresponding to the photons produced at the end of the crystal. On 
the other hand the antidiagonal width of (at (t)a(t')) is limited to maximum width of (a(t)a(t')) 
in the pump seed region, which is estimated as t p \J Lov /Ld- Therefore as soon as the crystal 
length exceeds Lov the spectrum of the SPDC stops shrinking. 



The Figures 2(e) and 2(f) show the CF in the later moments of evolution L = Lnl- All the 
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Figure 2. The contour plots of two Correlation Functions (CF) moduluses: \(a* (t)a(t'))\ and \(a(t)a(t'))\ presented on the 
t — t' plain. We draw a lossless case for the pump-SPDC overlap length Lqv = and dispersive length Ljj = 3Ljvl- 

The following pairs of plots present the evolution of the CF on the subsequent distances: L = Ljvl/20, ijv.t/3 and Ljvl- 





5 10 15 5 10 15 

t [r P ] t [r P ] 

(a) (aHt)a(t')) (b) (a(t)a(t>)) 



Figure 4. The contour plots of two Correlation Functions (CF) moduluses: | (a* (t)a(t')) | and \(a(t)a(t'))\ presented on the 
t — t! plai n. We dr aw the lossy case (the attenuation length La = 1/5Ljvl) corresponding to the other parametes as in the 
Fig, |2(c)] and [2(f)1 Lqv = Lnl/10, = 3Ljvl> L = ijvL- CF achieved a stationary state being unaffected during further 
propagation. 

above mentioned processes continue to shape the CF. The Fig. [3] presents contour plots of phases 
of (at (t)a(t')) and {a(t)a(t')) in the region of significant absolute values. They describe temporal 
chirp acquired by photon pairs in the propagation through the dispersive crystal. 

In the next step losses were added to the simulation. As we argued above for short propagation 
distances L <C Lov one can a Ppby ° ne mode solution given in Eq. (|2.4h . For longer distances 
two situations are possible. For high attenuation Ltvx > 2La and both of them much shorter 
than the overlap length Lov the CF will reach steady state values due to attenuation alone and 
dispersive effects will play a minor role. Therefore we show the results of the calculations only in 
the opposite regime, for La = 1/5Lat^ > Lov- I n this case the main role of losses is attenuation 
of the tail of CF which continuously unwinds from pump seed region. Thus the diagonal length 
of the CF is restricted to the length of the order of tpLa/ Lov as plotted in Fig. [2j The evolution 
at the distances L <ti La is very similar to lossless case. Then the CF stabilize until finally for 
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L 3> La they reach steady state when their tail is fully stretched and cut by attenuation. 



4. Evolution of the quadrature variances 



OPA produces light which exhibit squeezing of quadrature variances below the shot noise level. 
Such states demonstrate non-classical properties. Applying the Homodyne Detection (HD) to 
the light generated in SPDC we can directly observe the amount of noise in the mode of a local 
oscillator. Similarly to the previous sections, it will be instructive to start with the analysis of 
the single mode OPA. In a single mode HD we measure the quadrature x = A= (ae lip + a)e~ llp ^ , 
where <p denotes the phase difference between the pump and the local oscillator. In the SPDC 
light the mean value of x is zero. The variance of x which measures the quadrature noise equals 

(x 2 ) = 1 + (tfs) + Re{(aa)e 2 ^} (4.1) 

We are interested in the most squeezed and most antisqueezed quadratures. In the lossless case we 
can use (|2.4p and (|4,ip to retrieve a well-known result for the minimal and maximal variances in 
a single mode squeezing (x 2 ) = e ± 2L / L NL /2. The inclusion of losses characterized by attenuation 
length La changes this expression into: 

We see that depending on the relation between La and the nonlinear length Ljvl we deal with 
two different types of evolution. In the high losses regime La < Lnl/2 both minimal and 
maximal variances reach steady state. On the other hand in case of small losses La > Lnl/2 
the antisqueezed quadrature grows exponentially, while the minimal (x 2 ) decreases to a steady 
state value — the squeezing is limited. The above formula holds also in multimode case for 
crystal lengths much shorter than overlap length Lov- 

Let us follow with a multimode generalization of the previous discussion. We write down the 
expression for the quadrature x in the mode 4>lo(u) given by the macroscopic field of Local 
Oscillator (LO). 



x = —= I du> 

V2. 



^loM^M + 0Lo(w)a f (a;) 



(4.3) 



Variance of this quadrature (x ) can be obtained in a direct calculation: 



1 

2 + 



dtjdu/ 



4>LO (w) (a f (w)a (a/) ) 4>* LO <V ) 



+Re ((f) L o(u){a(uj)a(uj'))(p L o(uj')) 



(4.4) 



(^ 2 ) = 2 re f ers to the shot noise level. Let us note that first term under the integral does not 
depend on the phase of the (ftiOi while the second term containing (a(cj)a(u/)) is responsible for 
phase sensitivity of the noise. 

In most applications of SPDC we desire to select the most squeezed mode. Such a mode, 
determined by appropriate LO function 4>lo, minimizes (14. Ah for given CF (tf(w)a(u/)) and 
(a(uj)a(uj')) . Since we apply the matrix approximation of the correlation functions, the integral 
in (|4.4j) can be written as a quadratic form, multiplied both side by vectors representing real 
and imaginary part of (j)LO- The most squeezed mode is an eigenvector of the quadratic form 
corresponding to the smallest eigenvalue. This is described in details in the Appendix. 
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Figure 5. (a) The most squeezed function modulus (solid line) and phase (dashed line), (b) Histogram of the noise reductions 
versus mode number. 



We found the most squeezed mode function for the parameters used in the Sec. 3: Lqv 



1/10-Ljvx, = 2>Lnl, La = 1/5L/vl> L = L^l. The mode function presented in Fig. 5(a) 



is close to Gaussian. In the Fig. 5(b) we also present the squeezing versus mode number. The 
modes are found as subsequent eigenvectors of the quadratic form corresponding to (|4.4[) . which 
physically represent uncorrelated, squeezed quadratures of SPDC, as detailed in the Appendix. 
In a typical lossless case on the output of long crystals we obtain a number of modes exhibiting 
similar squeezing. Introduction of losses decreases squeezing most rapidly in higher order modes. 

Let us analyse a quadrature conjugate to x, = / du [— <p* LO (u>)a(ui) + (j>Lo{oj)w (uSj\ . 
In the lossless case x^/i is the most antisqueezed quadrature, and the product (x 2 )(x 2 ^ 2 ) = 1/4 
saturates the uncertainty bound. However in the presence of losses the product (x 2 }(x^ 2 ) grows 
during the evolution. In the lossy case this quantity can be treated as a measure of the impurity 
of the state and we will track its evolution. We have checked that in the lossy case the most 
antisqueezed quadrature differs from the quadrature conjugate to the most squeezed one. In 
addition to examining the evolution of (x 2 ) and (x 2 )(x^ 2 ) we will track the evolution of average 

number of photons in the SPDC (h) = f duj (w (u))a(u)) . 

Now we will be tracking the evolution of (x 2 ), (x 2 )(x 2 ^ 2 ) and (h) for different values of the 
parameters. Firstly, let us take the case with neglectable dispersion Lov^L^ — ► oo when the 
one mode solution (|2.4p can be applied. As it was discussed above, depending on the relation 
between La and Ljvl/2, different types of evolution occur. In the FigEKa) we present the big 
losses regime for La = 1/6 -Ljvl < 1/2 L^l. For L < La squeezing (x 2 ), number of photons 
(h) and the product of quadratures {x 2 )(x 2 , 2 ) grow. Afterwards losses stabilize all of them 
for L > La- Physically the attenuation balances a nonlinear amplification. Let us notice that 
regardless of the attenuation in the medium, one can always find squeezed quadrature in SPDC. 
Since losses continuously add a vacuum noise to the SPDC, and the mode becomes more classical 
as (x 2 )(x 2 ^ 2 ) grows, it is suprising that we always find some nonclassical properties in the output 
SPDC light. 

For La > Lnl/2 without any linear propagation, we encounter a stabilization of squeezing as 
in the previous case. Interestingly both (h) and (x 2 ){x?, 2 ) grow. The most squeezed quadrature 
remains constant but reduced below the shot noise level but its conjugate quadrature grows un- 
limitedly. When we deal with a more realistic situation where Lov is finite, the growth continues 
until L < Lov as presented in the Fig. [6^b). For longer crystals we observe a stationary state 
where the other two quantities (n) and (x 2 )(x 2 ^ 2 ) stabilize. 
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Figure 6. Plots presenting the squeezing, the average number of photons and the product of quadratures as the function 
of the propagation distance expressed in Ljvi units. In the plots we denoted the pump-SPDC overlap length Lqv an d the 
attenuation length La which values determine the type of the evolution, (a) The stationary state in the regime of big losses 
La < Lnl/2 with no dispersion Lov — * co (b) Stabilization of squeezing with the simultaneous growth of average number 
of photons and the quadrature product (x 2 )(x^^ 2 ) on the effective crystal length L < Lqv (c) Typical realistic case where 
La > Lqv leads to the stationary states of the correlation functions. Three stages of the evolution: nonlinear, linear and 
stationary, (d) Adiabatic attenuation of all tracked quantities from (c) with the pump dumping enabled. 



Typically the losses in the material are relatively small and we have La > Lov ■ This situation 
is presented in the Fig. [6|c). Here, the initial evolution of the SPDC light resembles the lossless 
case and it can be divided into two stages. During the first stage we observe the nonlinear 
growth of the average photon number and squeezing until L < Lov- I n the next stage squeezing 
stabilizes and the SPDC process generates many similarly squeezed modes. Thus the average 
number of photons increases linearly. The effects of losses can be seen in the third, stationary 
stage, when L > La- In this stage we reach the stabilization of all tracked quantities. In the two 
previously discussed cases, the stabilization emerged when the nonlinear interaction still played 
an important role. Here we deal with the attenuation of tails of the CF, described in details in 
Sec. 3 (Fig. [2]), where the multimode CF achieve the stationary state. 

In all of the discussed cases we neglected the second order dispersion Lp — ► oo. We have 
checked that the finite dispersion length Ljy does not influence the noise reduction in the most 
squeezed mode or the total average number of photons. Moreover, for a real media we usually 
have Lp^> Lov and the ratio between Lnl, Lov an d La plays dominant role. 

In the end we present what happens if we take into account both the pump and the SPDC 
light attenuation. Besides of the pump losses, in the Fig. [6](d) the other parameters are the same 
as in the Fig. Etc). Since the pump seeds the SPDC process its attenuation prevents reaching 
the balance between the losses and the nonlinear light generation. In the consequence the CF 
decrease adiabatically, together with the squeezing and the average number of photons. 
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5. Summary 



In this paper we studied the process of the multimode, one dimensional Spontaneous Parametric 
Down-Conversion in the lossy media. We utilized relatively simple approach to the problem 
consisting in analysis of the propagation of the first order Correlation Functions (CF) of the flu- 
orescence fields. These functions provide a complete information about the Gaussian multimode 
state of light. 

We wrote equations of propagation of the CF and expressed them using four characteristic 
length scales. We showed that the evolution of the CF can be qualitatively predicted, and the 
character of the solutions is fully determined by the ratios of the four length scales. In particular, 
in case of negligible dispersion the multimode problem simplifies to the one mode evolution, for 
which we found analytical solutions. We further differentiated the regimes of high and low losses 
for which we predict either a stabilization or an exponential growth of average number of photons 
in the SPDC light. The presence of losses stabilizes the most squeezed quadrature. However in 
the case of the small losses the antisqueezed quadrature grows boundlessly. 

In the multimode case we found that losses make the simple picture of separable modes of light 
not longer valid. However we are still able to single out uncorrelated squeezed quadratures. We 
described a simple procedure to find corresponding orthogonal modes. We tracked the evolution 
of the average number of photons in the SPDC. In addition we calculated the highest available 
squeezing and a product of the squeezing and the antisqueezing in the corresponding mode which 
we use as a measure of the impurity of the state. Examination of the development of these three 
quantities allowed us to divide the process into three general stages: an exponential growth and 
a squeezing due to a nonlinear interaction, a linear temporal expansion due to a group velocity 
mismatch and a stabilization due to losses. Those stages are well separated in a typical case when 
group velocity mismatch between pump and SPDC plays significant role on shorter distances 
than the attenuation does. 

We wish to emphasize that losses during the development of the spontaneous fluorescence light 
do not entirely destroy its quantum properties and the state always remains squeezed below the 
shot noise level. Under certain conditions this happens although both the number of photons in 
SPDC and the antisqueezing grow. 

Despite the fact the whole process contains the simultaneous effects of a squeezing, mode 
mixing and attenuation, we were able to divide the evolution into simple stages and predict it 
qualitatively. We find the description using multimode correlation functions very fruitful and 
promising. This approach could be also applied to other quantum problems where we observe 
simultaneous spontaneous scattering and losses. Possibilities include a nondegenerate three di- 
mensional OPA (fTl ). Raman scattering in atomic vapors 
BECs where spontaneously scattered atoms are observed 



l£\ and modelling of the collisions of 

3. 
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Appendix 



In the absence of losses the optimal mode minimizing (14.4p is the eigenmode of (a(uj)a(io')) with 
highest eigenvalue. Moreover we can find a set of uncorrelated modes squeezed independently. 
Unfortunately in the presence of losses neither is true. We have to minimise a general expression 
(j4.4|) by proper choice of the local oscillator mode 4>lo(^)- This is easily accomplished using 
linear-algebraic means. Firstly let us notice that the variance (i 2 ) in (|4.4p is real and can be 
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reexpressed in the following symmetrized way: 
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(x 2 ) = 1 + ff 



1 



&Lo{u){a\Lu)a(J))4>* LO {J)+ 



<& H(at(w)a(cy)>> L0 (u/) ) + Re (0 LO (^)(a(^)a(w , ))0Lo(^)) (1) 



As we apply a matrix approximation to the CF [a) (oj)a(uj')} , (a(uj)a(u>')} and vector approxi- 
mation to 0lo( w ) the integration in (pQ) becomes multiplication of a matrix by two vectors. The 
product is sensitive to the total phase of (pLoi^) an d therefore one has to represent it by a vector 
consisting of the real and imaginary part of (pLO'- 



Re(4>Lo) 
Im(0Lo) 



Then we can rewrite (pQ) in the following form: 

! ) = l/2 + $ T M$ 



ix 2 ' 



(2) 



(3) 



where M is the matrix of the quadratic form, expressed by the matrix approximations of CF 

(at (w)d(u/)) and (d(w)d(u/)): 



M 



Re ((a'a)) + Re ((aa)) , — Im ((aa)) + Im ((d^d)) 
-Im((dd)) - Im ((&*&)) , Re ((at a)) - Re((dd)) 



(4) 



Since Im ((d^d)) is antisymmetric, M is real and symmetric and its eigenvectors are orthog- 
onal. Therefore the correlation between two quadratures Xi and Xj corresponding to different 
eigenmodes and is zero: 



(xiXj) = *f M*, =0, i ^ j. 



(5) 
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